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ABSTRACT 



The advent of the missile and satellite era has brought 
with it the need for large and highly accurate radar tracking 
systems. Thus, the antenna tracking error generated by wind 
disturbances has become a significant consideration. This 
study presents a compensation technique which substantially 
reduces these disturbance errors without degradation of the 
system response to the tracking input. 

The study was proposed by Philco Corporation, Western 
Development Laboratories, Palo Alto, California. The proposal 
resulted from tests conducted on a production model of a 
satellite tracking system having a 60 foot, solid skin, high 
inertia antenna. A large portion of this thesis was accom- 
plished at the Philco Laboratories during the Summer of 1960. 

The authors wish to express their appreciation to 
Dr. M. P. Pastel of the U. S. Naval Postgraduate School for 
his assistance and counsel in the preparation of this thesis 
and to Mr. L. J. Harvey of Philco Corporation for his assistance 
in the development of the problem and for the generous facili- 
ties made available to the authors. 
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Introduction 



The purpose of this study is to present a compensation 
technique for reducing antenna tracking errors caused by 
wind-gust disturbances. The antenna system upon which this 
study is based is the position control system for a 60 foot 
radar antenna. This system meets the desired response to the 
tracking input. However, the tracking error caused by wind 
gusting is not acceptable. The compensation utilized must 
reduce this error without affecting the system response to 
the tracking input. This problem corresponds to having an 
antenna system installed at the tracking site which meets 
the designed input tracking response. Data showed excessive 
tracking error resulting from dynamic wind pressure on the 
solid skin of the large, high inertia dish. This error must 
be reduced by adding some type of compensation that will not 
modify the designed tracking input response. 

An error of a few seconds of arc can be excessive in 
determining the precise position of a high altitude orbiting 
satellite. To meet such stringent error regulating require- 
ments while, at the same time, achieving acceptable tracking 
input response using conventional compensation techniques 
poses a formidable problem. Compensating techniques require 
the sensing of the effects of the wind loading on the antenna 
system. Wind sensing might take the form of external devices 
such as sensitive annamometers . These might be placed in 
locations around the periphery of the antenna and at distances 
so as to determine in advance the load disturbance which the 
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antenna would see. Although the nature of a prediction 
sensing device lends itself to the development of the ideal 
regulation system, the accuracy of the prediction must of 
necessity be at least as accurate as the tolerance of the 
desired output response. To date, no wind sensing trans- 
ducers have been produced that measure accurately transient 
wind conditions over a sufficiently small increment of time. 

Even if such devices were feasible and available, their 
sensing would not necessarily reflect the actual loading 
of the antenna dish some distance away. A mountainous site 
might be expected to show large drafting and eddy-currents 
about the dish which would preclude the accuracy of any 
external sensing device to determine actual dish loading. 

A second and more reasonable approach might take the form 
of a sampling of the torque on the antenna drive shaft. 
Corrections would then be made in the control system based 
upon these samplings. This approach is complicated and 
expensive, and at least part time use of a digital computer 
would be required. 

A simpler and more direct approach would be to continuously 
compare the output position under external loading with the 
output obtained having no external load with the same tracking 
input. The difference of these outputs would then be a 
measure of the effects of transient wind loading on the 
antenna. Although this approach in no way predicts the loading 
effect, it does permit the accurate measurement of the effect 
which in turn may be compensated at a faster rate than that 
possible by the basic system. 
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The solution presented here will be developed in the 
above manner using as the basic approach the conditional 
feedback concept set forth in the American Institute of 
Electrical Engineers' publication, "Applications and 
Industry," July 1955, by G. Lang and J. M. Ham. The condi- 
tional feedback system permits design specifications on 
input-output response and on disturbance-output response to 
be met independently. 

To prevent repetition and to simplify the presentation, 
only the azimuth axis of rotation will be considered. Design 
criteria for the compensation will be developed based on 
selected arbitrary requirements for disturbance-output 
response improvement. A specific compensation will be designed 
using the developed criteria and the response will be 
determined for various disturbance inputs. A design technique 
is developed which facilitates the determination of this 
compensation . 

An investigation of the characteristics of a severe wind 
gust disturbance will be utilized to develop a disturbance 
input which will subject the compensated system to the most 
demanding conditions expected. 

An analog computer simulation will be made of the complete 
compensated system. Data will be displayed showing the 
effects of the selected compensation for various wind 
disturbance inputs. 
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2. Development of the basic approach. 

In the conventional single-loop feedback control system 
the input-output response, as well as the disturbance-output 
response, are dependent on the feedback. Thus, in most cases 
utilizing conventional feedback control, a compromise of the 
input response must be accepted if a satisfactory disturbance 
response is to be achieved. Fig. 2-1 shows such a conventional 
control system. 




Figure 2-1. Diagram of conventional control system. 



In Fig. 2-1, is the Laplace transform of some control system; 
G 2 represents the feedback compensation. R is the normal 
input signal to the system; U is the disturbance signal at 
the output of the system; C is the final system output. 

The transfer functions associated with Fig. 2-1 are readily 
developed as follows: 



C = Gj^CR - G^C) + U 



C(1 + G^G^) 
C 



Gj^R + U 

Gi ^ 

1 * GjG2 * 1 * 



(2-1) 



Thus, the feedback G 2 affects both input and disturbance 
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responses. However, the problem treated in this thesis is 
one in which the effects of disturbance must be compensated 
without compromise of the given input-output response. A 
compensation technique meeting these requirements has been 
presented in ’’Conditional Feedback Systems" by G. Lang and 
J. M. Ham [l] and will be used as the basis for the approach 
of this thesis problem. 

A basic configuration for a linear conditional feedback 
system is shown in Fig. 2-2, where , G 2 , R, U, C, are defined 
as in Fig. 2-1. 




Figure 2-2. Basic conditional feedback control system. 

The output C of Fig. 2-2 is 

G.(l + BG ) 1 

^ " 1 + Gj^G2H ^ 1 + G^G2H (2-2) 

Let B in equation (2-2) be defined as 



B = Gj^H. 

Equation (2-2) then becomes 

C = G.R + U. 

1 . G^G2H 



(2-3) 



(2-4) 



5 







I 



Letting U and R be alternately zero, 




G 



1 



(2-5) 



and 




1 



(2-6) 



R=0 



1 + 



Equation (2-5) shows that, when B is defined by equation 
(2-3), the input response is dependent only on Thus, G 2 H 

can be utilized to modify the disturbance response without 
affecting the input response. The output sensing device 
would normally be described by H. If it is assumed that 
H = 1, then B becomes equal to Gj^ (equation (2-3)). The 
purpose of the function B is shown by considering the system 
operation in Fig. 2-2. The feedback signal (b) is equal to (a) 
since H = 1. With B defined as equal to G^^ , (d) is equal to 
(a) with U equal to zero. Thus, the difference signal (e) is 
equal to zero and no compensation signal is developed. 

Clearly, (d) is equivalent to the desired output response and 
thus B is called a "reference model." The reference model 
is an undisturbed representation of the transfer characteris- 
tics of the antenna drive system. This conditional feedback 
approach permits the use of the feedback function G 2 to 
accomplish the desired disturbance compensation without 
change to the original input response. 
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3. Development of basic system. 



The basic system investigated in this thesis represents 
the azimuth axis position control system of an actual land- 
based tracking radar having a sixty foot diameter antenna. 

The system consists of a hydraulic drive with tachometer 
feedback and suitable reduction gearing, together with cascaded 
electronic compensation and position feedback loop. 

Prior to any possible investigation of the basic system, 
the hydraulic drive system transfer functions had to be derived. 
This hydraulic system consisted of a servo control value and 
amplifier, a D. C. motor-driven constant displacement pump and 
associated hydraulic drive motor, and reduction gearing. The 
hydraulic system was enclosed with a tachometer feedback loop. 

The block diagram of this hydraulic drive system is shown 
in Fig. 3-1. 



error signal 




Figure 3-1. Hydraulic drive system. 

The transfer function is defined as the joint function 

of the servo amplifier, servo control valve, and hydraulic 
drive motor. Neglecting the tachometer feedback (G^), the 
output position angle to the reduction gearing is (Oj^), where 

©1 = Gj^X (3-1) 

and X is the compensated error signal of the basic system. 
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The block diagram of the hydraulic drive system with normal 



antenna loading and external load torque (T^) due to wind is 



shown in Pig. 3-2. 
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Figure 3-2. Hydraulic system with synthesized loading. 



Letting = 0, 

T = T + T 
1 m _a 

N 



T 

m 




(3-2) 



where T is the torque due to hydraulic motor inertia (J ) 
and is defined as 

T„ = J s^0. . (3-3) 

m m 1 

The load torque (Tj^), as shown in Pig. 3-2, is 

A = * A = • (3-4) 

©1=9^, (3-5) 

TT 



where N is the gear reduction ratio of the hydraulic system. 
Substituting equation (3-5) into (3-4) gives 

A = A'A - ®L> = <A=^ * A">®l " A ■ (3-«') 

N 
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or 



(3-7) 



K, 



®L = 



N(Jj^s^ + fj^s + Kj^) 



8l - 



(Jl= " " ''l) 



Substituting equation (3-7) into (3-6) gives 

.2 



K 



K. 



- 









(JlS + f^s + K^) 



(3-8) 



The flow of oil out of the control valve (Q^) is considered 
to be divided into three parts, that lost to system leakage 
(Q, ), that lost because of compressibility of the oil (Q ), 
and that actually used to drive the hydraulic motor (Qj^)* 

Qv = Ql » Qo * Qm «-9) 

The flow out of the control valve is a function of the input 
error voltage (X) within some constant (K^) and is defined as 
Q = K X . (3-10) 

^v V 

The constant (K^) is defined to include the control valve 
constant as well as the servo amplifier gain constant. Brown 
and Campbell in "Principals of Servomechanisms" [2j define 
the separate oil flows of the system as follows: 

Ql = L P . (3-11) 

where L is the leakage coefficient and P is the pressure 
differential of the system. 

Qc = I sP , (3-12) 

where V is the volume of oil under compression and B is the 



bulk modulus of the oil. 
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(3-13) 



Q 



m 



dm , 



where dm is the motor displacement and 0^ is shown in 

Pig. 3-2. Dimensions of all constants are listed in the table 

of constants on page 12. 

Making substitutions of these flows into equation (3-9) gives 
K X = dm s0, + (L + ^s) P. (3-14) 

V Id 

The pressure differential (P) of the system is related to the 
torque output of the motor (T^) as follows: 

= dm P . (3-15) 



Substitution of equations 
(3-2) gives 

T, = dm P = J s^0, + 
1 ml 



(3-15), (3-3), (3-4), 



(JlS^ + f^s + Kj. ) 



N(Jj^s^ + f^s + Kj^) ^ 




into equation 



+ 1 



0 



1 



(3-16) 



or 



J s^ K, 

P = — — 0^ + — - 



(JlS + fj^s) 



dm 



dm (Jj^s^ + fj^s Kj^) ^ 



K, 



dm N ('IlS + fj^s + Kj^) 



(3-17) 



Substitution of equation (3-14) gives 
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1 




K X = dm sO, + (L + Vs) 
'' 1 B 



j e, K, (j.s^ + f, s) e, 

m 1 ^ L L L 1 

dm 



dm N^CJj^s^ + fj^s + Kj^) 



(L * |s) 



dm N (Jj^s^ + fj^s + Kj^) ^ 



(3-18) 

If the external torque (T^) is assumed equal to zero, then, 

V 

solving for ^ in equation (3-18) gives 



®1 ■ ^ 



0 . 



dm s + (L + ^s) + (L + |s) — 2 



(JlS + f^s) 



dm 



dm N (Jj^s +fj^s+Kj^) 



Since G, = 1 , 



G,. = — = 



K 



‘dm + (L + Vs)J s + (L ♦ Vs)K, 

B -S!- B — 



‘V * 



dm 



dm + Kj^) 



(3-19) 

Since the system stiffness factor (Kj^) is relatively large, 
equation (3-19) reduces to 






BK^ dn / V(J^ 



s + (13 + 
V 



) s . LBf^ 



dm^ B 



N^(J + J, ) N^(J +J, ) V(J + J, ) 
m L m L m L 



n2 



n2 



n2 _ 
(3-20) 

The friction coefficient of the antenna drive is considered 
very small with respect to the load inertia constant, thus. 



N"'(J + J, ) 

m L 



= 0 , 



N" 
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and equation (3-20) reduces to 



B dm / V(J, + N^J ) 
V ' L m 



S = 



+ LBs + dm^ B 



V 



V(J, + N^J ) 

L m —I 



(3-21) 



The hydraulic constants associated with this system 
are as follows: 



N 

K 

J 

J 



V 

L 



B 

L 

V 



dm 

Substi tuting 



= 2160 
= 0.895 

6 2 

= 4.2x10° Ib-in-sec (referred to antenna) 

2 

= 0.09 Ib-in-sec (referred to motor) 

= 0.27x10^ psi 
= 60.5x10"'^ in^/sec/psi 
= 3.9 in^ 

= 0.215 in /radian 

these constants into equation (3-21) gives 



S = 



1.345x10 



s(s^ + 42s + 3228) 



(3-22) 



Equation (3-22) is the transfer function of the 
uncompensated hydraulic system and was used in the initial 
development of the complete basic system which follows. The 
transfer functions of the complete basic system chosen initially 
for investigation are shown in Fig. 3-3. The values shown for 
the individual transfer functions were suggested by the 
manufacturer as representing the azimuth control system for 
the antenna under investigation. The open-loop transfer function 



(F^) and the closed-loop function (F^) for the system of 
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= Cascade compensation 
Gj^ = Hydraulic drive 

G£ = Tachometer feedback (tachometer constant = .086) 
N = Gear ratio (2160) 



Figure 3-3 

Block diagram of initial basic system. 
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Fig. 3-3 are as follows: 



p = 14.96 X lo'^Cs 1.7)(s 20) 

° s^(s + 10)(s^ + 42s + 4038) 

14.96 -f lO'^Cs 1.7)(s + 20) 

^ s^ + 52s"^ + 4458s^ + 1.8988xl0^s^ + 3.24xl0^s 



(3-23) 



5.08x10^ 

(3-24) 
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4. Analysis of basic system. 



In order to ascertain the stability of the system of 
Fig. 3-3, a Routh's criterion test was made of the closed- 



loop transfer function 

Routh's criterion test: 

1 4.458x10^ 3.24x10^ 

52 1.8988x10^ 5.08x10^ 

810 3.14x10^ 0 

-1.172x10^ 5.08x10^ 0 

3.49x10^ 0 0 

5.08x10^ 0 0 



The system as given was believed stable, but Routh's 
criterion test did not verify this assumption. Accordingly, 
a Nyquist plot was drawn of the system transfer function 
in order to observe the degree of instability and is shown 
in Fig. 4-1. Within the accuracy of the plot, the curve 
appeared to pass through the Nyquist point indicating marginal 
stability. For further investigation of the system stability, 
a Bode plot, using asymptotes only, of F^ was constructed and 
is shown in Fig. 4-2. Although the plot crossed the zero db 
line with a negative six db slope, the complex poles of F^, 
having a natural frequency of 63.5 rad/sec, caused a hump 

in the plot which, if the magnitude were sufficiently large, 
would recross the zero db line and result in system instability. 

In order to estimate the maximum amplitude of the hump, 
that part of the system transfer function (^) resulting from 
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the complex poles can be written in vector form as 



F = 



= Me 



jN 



(4-1) 



(1 - j2S^ 



w 



w, 



w. 



where: w^ = break frequency of the complex poles 

S = damping coefficient 
M = IfI 

A 

N = phase angle of F 

The maximum absolute value of ^ (M ) is the maximum 

max 

amplitude of the hump. Letting ili = u and substituting into 

'''o 

equation (4-1) gives 



( 1 - u^) + j2 S u 

or 

M = ^ . (4-2) 

\|(1 - u^) + 4 S ^u^ 

Differentiating equation (4-2) with respect to u gives 
dM = -2(1 - u^)(-2u) + 8 § ^u 

2 [(1 - u2)2 . 4^ VJ i [(1 - u2)2 . 48 V] 

(4-3) 

Setting equation (4-3) equal to zero and solving for u gives 
u = (1 - 2 S ^)'^ . (4-4) 

Substituting equation (4-4) into (4-2) yields 

^Chestnut and Mayer, Servomechanisms and Regulating 
Systems Design, Vol I, pp 312, September 1953. 
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1 



(4-5) 



2 « ( 1 - » 

Equation (4-5) is an approximation of since it was 

developed on the basis that the frequency at was equal 

to the break frequency (w^). However, these two frequencies 
differ slightly. 

The damping factor (S ) can be computed from equation 
(3-23) as 

2S = 42 , 

where w = (4038 )t , or « = 0.33. 
n 

Therefore, from equation (4-5) 

M = 1.7 = 4.62 db . 
max 

The maximum amplitude of the hump, 4.62 db , was added to the 
amplitude at the complex pole break frequency in Fig. 4-2 and 
resulted in a system amplitude at this frequency of 0.12 db. 

Although the above calculation is not exact, it at least 
shows the hump quite close to the zero db line and thus, the 
system has a very small stability margin at best. A decrease 
in system gain would improve the stability by lowering the 
hump below the zero db line, but as only 5.5 db can be dropped 
before the slope crossing the zero db line decreases to -12 db, 
this approach is not satisfactory. 

The above marginal stability was a result of the location 
of the complex poles of . These poles in turn were placed 
by the location of the complex poles of the hydraulic drive 
transfer function (G^) and the tachometer feedback (G£). Near 
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cancellation of these latter poles was achieved by insertion 

of a complex-zero compensator (G^) in series with G^, 

alleviating the instability problem. The design of the complex- 

zero compensator appears in Appendix I. 

Equation (3-22) has shown the complex poles of to be 

(s + 42s + 3228). The form of the transfer function of the 

complex-zero compensator (G^) is shown by equation (1-9) in 

Appendix I. Thus, the G^ necessary to cancel the complex poles 

of G, is 
h 

G = + 42s + 3228). 

o 322o 

For the purpose of this study, complete cancellation of the 

complex poles of G^ was assumed. In the actual application, 

some residue would remain since equation (1-9) does not 

exactly represent the transfer function of the compensator 

in Fig. I-l of the Appendix. 

Placing G in series with G. inside the tachometer 
° o h 

feedback loop results in a new open-loop transfer function, 

F = 37 (s^ 1.7KS . 20) ^ 

° (s + 10) 

Fig. 4-3 shows the Bode diagram of the function of 
equation (4-6). A comparison of this Bode with Fig. 4-2 shows 
that the desired cancellation has been achieved without change 
of the system gain or location of the remaining break points. 
Thus, the system response remains the same with the exception 
of the undesirable effect of the complex poles. 

At this point, a limitation was placed on the system 
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cascade compensation (G^) . In order that components of less 

stringent design might be utilized in , the cascade gain of 

4 r 3 

2.405x10 was limited to a value of 1x10 . In order that the 

system gain of equation (4-6) be retained, servo amplifier 

gain of the hydraulic system as well as tachometer feedback 

gain were modified through the following development. 

The open-loop transfer function of the basic system of 

Fig. 4-4 is 



F 

o 




G G. 
o h 



1 + 



G G. G. 
o h f 



(4-7) 



Equating the system gains of equation (4-7) with the desired 
system gain of equation (4-6) gives 




K K. 
o h 



1 + 



K K, 
o h t 



37 , 



(4-8) 



where K represents the gain constants of 
transfer functions of equation (4-7) . 
in equation (4-8) gives 

Vh _ 37x2160 _ on 

1 + K K, K. 1000 • 

o h f 



the subscripted 
Substituting = 1000 

(4-9) 



Substituting = 0.7x0.086, where 0.086 is the tachometer 
constant, into equation (4-9) gives 

(1 - .7x.086)K K. = 80 
o h 

or 

= -20.95 . (4-10) 
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Since all transfer function gains must be positive, equation 
(4-10) indicates that the tachometer gain must be changed in 
order to realize the desired overall system gain. In order 
that be a positive number, equation (4-10) demands 

(1 - 80K^) i 0 . (4-11) 

Since the tachometer constant (.086) can not be changed, the 

tachometer gain must be changed so that its value is less 

that 0.1455 in order that K K, be positive. To ensure this 

oh 

positive value, a tachometer gain of 0.135 was selected. 

= .135X.086 

and 

(1 - .135x.086x80)K^Kj^ = 80 , (4-12) 

or 

K K. = 1150 
o h 

where K = j;- ; therefore 

O ^ ^ o 

= 1150x3228 = 3.72x10^ . (4-13) 

Equation (3-21) defined as 

B dm K 

S ^ V(Jj^ + |^2j^) “ 1.528x10 . (4-14) 

Substitution of equation (4-13) into (4-14) gives 
\ = 245 

where was previously defined as including the control 
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valve constant as well as the variable servo amplifier gain. 

The transfer functions of the basic system reflecting 
the above changes are 



g ^ 3.72x10^ 

^ s(s^ + 42s + 3228) 



G = (s^ + 42s + 3228) 

o Jz2o 



= ( .135x.086)s 

„ 1000 (s + 1.7)(s + 20) 

^c s(s + lO) 

The new open-loop transfer function (F^) was shown in 
equation (4-6). The new closed-loop transfer function was 
derived : 

P ^ 37 (s + 1.7)(s + 20) 

^ (s + 1.74)(s^ + 45.26s + 724.4) 



This is the basic system transfer function referred to here- 
after as G^ and is shown in Fig. 4-4. 

As was shown in the development of the basic approach, 
the conditional feedback compensation effects only the 
disturbance response of the basic system. Fig. 4-5 shows 
the response of the basic system to a unit step disturbance 
input (U) , where, 



The response curve of Fig. 4-5 was plotted using the transient 
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Figure 4-4. Block diagram of basic system with 
complex zero compensator and gain modifications 
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response equation determined by the graphical residue 
technique. The details of this procedure are included in 
Appendix II. The maximum overshoot in Fig. 4-5 is 10.3% 
and occurs at t = 0.10 seconds. Fall-time, defined as the 
time interval between ten percent and 90 percent decay during 
the initial cycle of response, is 0.037 seconds. The condi- 
tional feedback compensation to be attempted will try to reduce 
significantly this fall time. 

Fig. 4-6 shows the frequency response of the basic 
system to a disturbance input of from 0.5 cps. to 40 cps . 

It is noted that as the frequency of the disturbance input 
signal is increased, the regulating effect of the basic 
system deteriorates. In order to compare this frequency 
response with that resulting after conditional feedback 
compensation, the system bandwidth is defined as that 
frequency above which the input disturbance is attinuated 
by less than 10.62 db, or where M is equal to (1 - .707). 

In Fig. 4-6 this value of amplitude is reached at a frequency 
of 2.10 cycles per second. The corresponding phase angle 
is (-)78 degrees. 
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5 . Development of compensation technique. 

The following is a derivation of the conditional 

compensated disturbance-output transfer function 

together with a presentation of a compensation technique. 

Fig. 5-1 is a block diagram showing the inputs and transfer 

blocks of the proposed conditional feedback compensated 

system, where 
F 

G. = . closed loop transfer function of the 

^o basic system 

F^ = the open loop transfer function of the basic system 
G 2 = compensator to be determined 
R = tracking input information 
U = load disturbance information 
C = system output 




Figure 5-1. Block diagram of proposed conditional 
feedback system. 
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The transfer function of the output (C) with respect 
to R and U is developed from Fig. 5-1. The relations below 



follow from Fig. 5-1. 

C=U+c=U+(a- C)F^ 

c = bF = (a - OF 
o o 

b = a - C 
d = RGjl - C 
e = RG 2 Gj^ - CG 2 
a = R+ e = R+ RG 2 Gj^ - ^^2 

Combining these relationships gives 

C = U + (R + RG^G. - CG„ - OF . 

Z 1 Z o 

Rearrangement of equation (5-1) gives 



C 



U 

1 + F + F G„ 
o o 2 



+ 



(1 + GiG,)F 
L Z O 

1 + F + F G^ 
o o 2 



F 

Since 1 + F = 

o (ji 



then equation (5-2) becomes 



C = 



U 



+ (1 Gi*^2)Pq 



^ ^ "oS 



7 =^ + F G 



o 2 



C = 1 u ^ (1 H- GiG 2 )F^ 

^(1 + Gj^G2) ^(1 + Gj^G2) 

^ “ (1 + F )(1 + G.G,) ^ ^1^ 

O 1 z 



(5-1) 



(5-2) 



(5-3) 
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Therefore, the output due to load disturbance (U) is: 



[ 



1 



1 

+ 



P 

o J 



1 J 



U. 



(5-4) 



The output due to the tracking input (R) is: 



Gj^R . (5-5) 

The input-output transfer function is the same as that 
for the uncompensated basic system. Thus, this input-output 
response remains unchanged with the addition of the compensation. 

Root-locus techniques are to be used to select a G 2 which 
will give a desired disturbance-output response. The disturbance- 
output transfer function may be written in the following form 
to facilitate the application of the root-locus technique. 



t * ^o] t * '=1' 



O 2 . 



2J 



1 + P + P G^ 
o o 2 



+ N N_ + N 
0 2 o 2 o 2 



(5-6) 



where : 




P 

o 

1 + P 

o 




N 
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The denominator of the above function requires two root loci 
in order to determine the final closed-loop roots of the 
disturbance-output transfer function. Since it is desired 
to investigate the effects of several values of the compensator 
it seemed advantageous if a function, closely representing 
the disturbance-output function, could be utilized which 
would require only one root-locus plot to place the roots of 
the closed-loop function. An examination of the function 



1 + 



^ 1^2 



shows it to be such a function. 



This function can 



be rewritten as 



1 + P( 



^ 1^2 



1 + 






1 + F 



F G., 
o 2 



^ 1^2 



^ 1^2 



where 




(5-7) 



The roots of the function of equation (5-7), hereafter referred 
to as the "representing function,” are the same as that for 
the desired disturbance-output function of equation (5-6). 
Furthermore, this function is of a form which requires the 
use of only one root locus in order to place the roots of the 
desired closed-loop function. Thus, if a direct correlation 
between the representing function and the desired function can 
be shown, the problem of selecting G 2 through simple root- 
locus techniques is made relatively easy. The development 
of this correlation follows. 
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The representing function of equation (5-7) can also be 



written as 



D D_ + N 
o 2 o 2 



1 . G,G2 



D + N -*■ N 
o 2 o 2 o 2 



For the basic type 2 system given, takes the form: 

2 

s ( r + 1)( r + 1) 

and represents the denominator of the open-loop transfer 
function of the basic system. 

Assuming the compensator (G 2 ) takes the form: 

s'^Cr ,s + l)(r s + 1) 

i ^ , m > 0, 

(r 3S + i)(r ^s + 1) 

then the numerator of the representing function, (D D_ + 

N^D 2 ) , is some polynomial in s of the form 

(As” + Bs”'^ +....+ D) . 

The denominator of the representing function is also of the 
same general form, (A's^ + B's^”^ +....+ D*). Since 
no factored s term exists in either the numerator or denominator, 
the response to a step input of this representing function will 
include a residue at the "input” pole inserted at the origin 
and, therefore, will have a steady-state error constant. 

The desired disturbance-output function has been shown 
in equation (5-6) to be: 

^ 0^2 

D N D., N ~ N ~ ’ 

o 2 o 2 o 2 
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s^(r + i)( r 2S + 1) . . . . (r 3S + i)(r + 1) 



A’sP + B'P“^ +....+ D’ 



Since this form does not produce a residue at the origin 
of the s plane due to a step input, the time response has 
no error constant. 

The time response of the representing function to a 
step input is of the form: 



approximation. As time increases, the value of the output 
of the representing function will tend to be larger by an 
amount thaji that expected for the desired function in which 
is equal to zero. 

If there exists no phase difference ( T ) between the 
representing function and the desired function, then to a 
second order approximation, the correlation would simply be: 



This, of course, is not the case. Some phase difference 
between the two functions must exist since they are not 
identical. It is this phase difference ( (T ) between the 
functions which in addition to the constant error term, 
decides the similarity between the two functions. Some 
correlation with respect to this phase difference is necessary. 

An inspection of the desired transfer function. 



Cl + €30 *^1^ sin(w^t + f ) , to a second order 
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shows that the zeros of this function are made up of the 



poles of F and G„ . An inspection of the representing function, 
o ^ 



1 + Gj_G2 



) = 



D»(N^ + D^) 
2 o o 



D + N + N ’ 
o 2 o 2 o 2 



shows this function to have zeros that are made up of the 
poles of Gj^ and G^. Both functions contain zeros which are 
the poles of G 2 . In addition, both functions have the same 
roots. The only difference between the time response of these 
functions must be due to the difference between the placement 
of the poles of F and the poles of G, , (N + D ), the 
remaining zeros of the two functions. To a second order 
approximation, the phase difference (T) can be determined by 
noting the placement of the poles of F^ and G^ with respect 
to the complex pole of either the representing or desired 
function, since they are identical. The phase difference ( T ), 
determined using graphical residue methods, is as follows; 

y = Sum of the angles from the poles of F to 
the complex pole of the representing ° 
function minus the sum of the angles from 
the poles of G^^ to the complex pole of 
the representing function. 

(5-8) 

Since the response of the representing function to a 
unit step input, to a second order approximation, is of 
the form 

Cf + sinCw^t + y) , (5-9) 

then the response of the desired function to a unit step 
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input must be, 



C^e" ^^sin(w^t + V + 7T ) . (5-10) 

The development thus far permits the determination of 
the correlating parameters, phase difference ( 'S' ) and steady- 
state error constant (Cj^) . Techniques for utilizing these 
parameters to predict the desired response from the representing 
response are presented in the following development. 

A comparison of equations (5-9) and (5-10) show that 
they differ by a constant, a proportionality factor, and a 
phase difference (S'). The value of can be determined to a 
second order approximation and is shown below. 




The representing function, referred to hereafter in this 
section as R.F., has pole and zero locations as shown in 
Fig. 5-2. Applying graphical residue techniques to Fig. 5-2 
for a unit step input gives 
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= Distance to zeros of R.P. from the origin 
divided by the distance to the complex poles 
of R.F. from the origin. 



Since is chosen such that R.F. approximates a second order 



system, the zeros of R.F. are equal to the poles of Gj^, thus 

2 

(5-11) 



r = 
^1 



Wn(Gi> 



where is the distance from the complex poles of G^^ to 

the origin and W^(R.F.) is the distance from the complex poles 
of R.F. to the origin. Equation (5-11) shows, for example, that 
if the compensated response is to be three times as fast as 
that of Gj^, then W^(R.P.) should be approximately three times 



that of this example. 



^ (i)^ = .11 

Thus, for a significant improvement in response time, will 
be small. If the phase difference (If) is also small, C 2 is 
approximately equal to since both equations (5-9) and 
(5-10) must equal unity at zero-plus time. With both and T 
small the response to a unit step input of the desired function 
will be similar to that of the R.F. For small and T the 
overshoot of the desired function may be predicted. The 
overshoot of the R.F. to a second order approximation can be 
determined from the root locus using simple second order 
techniques. The value of is then added to this estimate 
of R.P. overshoot and is used as the prediction of desired 
function overshoot. 
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Equation (5-10) shows the phase angle of the desired 
function to be equal to the sum of V and Y . The value of y 
together with its sign was computed using equation (5-8). 

If a prediction of the desired response based on the R. P. 
response is to be made, some knowledge of S' is necessary 
in order that the effect of ^ , the phase difference, may be 
known. For example, let S' be a second quadrant angle and r be 
positive. During the initial period of response, the sine 
term of equation (5-10) will take on smaller values than 
will the sine term of equation (5-9). This is to say that 
for small C^, the response of the actual desired function 
will be somewhat faster than that of the R. P. Similar 
examples can be given for negative y or for y in different 
quadrants. Suffice it to say that a knowledge of the quadrant 
location of y is necessary. The shaded area of Fig. 5-2 defines 
an area where the location of the complex poles of R. P. 
produce a v always in the second quadrant. This area is 
defined in terms of known parameters. It is that area which 
lies to the left of the complex pole of where any complex 

root of R. F. will have a damping factor ) greater than or 

equal to that of G^, and will have a damped frequency of 
oscillation (w^) greater than that of the complex poles of 
G^. This area, in general, represents the possible locations 
of the complex pole of R. F. which will produce a faster 
response than G^^ but will not effect an increase in over- 
shoot since the damping factor is at least as great as that 
for G^ . Thus, the most probable location for placement of the 
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poles of R.F. lies in an area in which y is known to be a 

second quadrant angle. For complex poles of R.F. in this area, 
the determination of a positive y will predict a faster 
response and a negative y will predict a slower response 
of the desired function as compared to the representing 
function (R.F.). As the absolute value of r increases, the 
responses of the two functions become more dissimilar. 

Accuracy of the prediction of the desired response 
from the representing response is gained with experience in 
using this technique. The steps utilizing the "representing 
function technique" for determination of the conditional 
feedback compensation G 2 are summarized below. 

1. Specify the desired disturbance response of the 
compensated system with respect to the original 
basic system response. 

2. Utilize the root-locus method to determine a G 2 that 
will place the closed-loop poles of R.F. in a position 
such that desired specifications are approximately met 
and the R.F. response is approximately of second order. 

a. Place poles in that defined area in which 

y is a second quadrant angle. 

b. Choose of R.F. complex poles such that 

Wj^CR.F.) 

^ - Q — ^ — is equal to the desired improvement 
n 1 ' 

in the response time. 

c. Pick a damping factor to meet desired overshoot 
requirement . 
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3. 



4. 



5. 



Determine correlation parameters. 

a. Compute = 

b. Compute Y using rule of equation (5-8). 

Make first prediction of desired function overshoot. 

a. Using second order approximation technique, 
estimate overshoot of R.P. 

b. Sum the estimated R.P. overshoot and C^^. 

Predict the desired function fall-time. 

a. Estimate fall-time of R.P. response as, 



Wn(Gi) 

W^(R.P.) 



W (G ) 

fall-time (R.P.) = fall-time (G^^) x ( R . p ) 

n 

(5-12) 



b. Predicted fall-time of the desired function 

response is the above estimate modified by r . 
Positive Y predicts less fall-time and negative r 
predicts greater fall-time than that of the R.P. 
response. The degree of modification is propor- 
tional to the absolute value of V . 

6. Make second prediction of desired function overshoot. 

a. Increase first predicted value if T is positive 
and decrease if negative. The degree of modifica 
tion is proportional to the absolute value of "JT . 

7. Compare predicted response with specified desired 
disturbance response. If satisfactory, proceed to 
verify response of the actual system including the 
G 2 conditional feedback compensator. 
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6. Determination of conditional feedback compensator. 

The representing function, ^ ^ , will be used to 

1 ^ 1^2 

determine the specific compensation through root-locus 
techniques in the manner outlined in Section 5. Observation 



of the response of the basic system to unit step disturbance, 
Fig. 6-1, shows it to have the following characteristics: 
Fall-.time: .0375 sec. 

Overshoot: 10.3% 

Steady-state error: 0 

Number of oscillations: 1 

In order to effect a significant reduction in the track- 
ing error of the system due to transient disturbances, the 
transient characteristics of the representing function were 
arbitrarily chosen as follows: 

Fall- time: Decrease by 300%> to 5007o. 

Overshoot: No increase over that of 

basic system. 



Steady-state error: 



Small , 



Oscillations: Essentially the same as 

for basic system. 

The procedure followed in using the root-locus method 
was to choose a G.. which would effect an increase of W (G, ), 
the natural frequency of the basic system, by an amount of 
300% to 5007o. At the same time, the real roots of the basic 
system were compensated in such a way as to minimize the 
residues of the real roots of the resulting representing 
function (R.F.). Having done this, the R.F. approximated a 
second order system and its time response was generally known. 
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^ I 







I 

I 



I 




The simple correlations developed in Section 5 would then 
predict the desired function response. 

The basic system closed-loop transfer function (Gj^), was 
given in Section 3 as, 



P _ 37(s + 1.7)(s + 20) 

1 Cs + 1.74)Cs + 22.63 + jl4.53)Cs + 22.63 - J14.53) 

therefore , 

^ (s+1.74)(s+22.63+ jl4.55)(s+22.63- J14.55) 

l+Gj^G^ ” (s+1.74)(s+22.63+ jl4.55)(s+22.63- J14.55) + ' 



where G^ 




+ 37(s+1.7)(s+20) N ’ 

^ 6 - 1 ) 



Several values of G 2 were utilized in an attempt to 
meet the time response requirements arbitrarily chosen for 
R.F. The following are representative of the compensations 
attempted . 



K 

( s+a) 

K(s-i-b) 

(s+a) 

K(s+b) 
( s+a) 



a lag function, 
b > a, a lag function, 

a>b, a lead function. 



Ks , a derivative function. 

The G 2 compensation which gave the desired position of 

the roots of the representing function (R.F.) was, 

^ _ 21.06(s+96)(s+200) 

2 (s+20)(s+1000) ’ 



a lag-lead compensating function. Fig. 6-2 shows the root- 
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^4 



locus of the complex roots of the R.P. utilizing this ^ 2 ’ 

The complex roots are shovm at 

s^ = -80-J52 ; S 2 = -80+J52. 

The real roots of R.F. were solved numerically and were found 
to be : 

s^ = -1.712 ; s^ = -20 ; s^ = -1660 . 

The natural frequency of the basic system, is 26.9 

rad. per sec. The natural frequency for R.F., W^(R.F.), is 

95.35 rad. per sec. This is an increase of 350% over '^^^^1^’ 

Section 5 has shown the form of the representing 

— 80 1 

function to be + ^ 2 ^ sin(52t + T) 

for the complex roots of R.F. determined above. The phase 

angle ( ) lies in the second quadrant since the complex poles 

are located in the area that insures a phase angle in the 

second quadrant. The steady-state error constant (Cj^) was 

approximated using equation (5-11) to be .079. Since the 

roots of this function are the same as those of the desired 

function, the latter function, to a second order approximation, 

- 80 1 " 

takes the form, sin(52t + V + ■»■ ) . 

In order to apply equation (5-8) for determining the 
phase difference ( V), the poles of F^, G^, and the complex 
poles of R.F. were plotted in Fig. 6-3 and T was measured 
as +13 degrees. 

The R.P., to a second order approximation, should 
give a fall-time of about 350% faster than that of the basic 
system, or applying equation (5-12), 

Fall-time (R.F.) = 26 . 9 rad/ sec x 0.0375 sec = .01 sec. 

95.35 rad/sec 
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Since Y is positive and a relatively small angle, it. is 

expected that the fall-time of the desired function response 
will be slightly less than that estimated for the R.F., that 
is, slightly less than .01 sec. 

The overshoot of the R.F., to a second order approximation, 
was 0.008 . This value added to the of 0.079 gives a first 
prediction of desired function overshoot equal to 8.7% . 

Since Y was found to be positive, the overshoot of the desired 
response should be increased to some value slightly greater than 
that of the first prediction of desired function overshoot, that 
is, slightly greater than 8.7% . Thus, the predicted response 
of the desired function falls within the parameters selected 
for the new compensated system. 

The desired transfer function is. 



(1+F^)(1+G.G,) 
o 1 z 



D 

o 

(D +N ) 



D„(D^ + 
z o 






(D D, + 
o 2 



NDT^TlOirT 
o 2 o 2 



The roots, (D + N D., + N N.,), have been found to be, 
o 2 o 2 o 2 

( s+l. 712) (s+20)(s+ 1660) (s+80+j52)(s+80-j52) . 



Therefore, the desired transfer function is, 

s^(s+l .74)(s+10)(s+20)(s+1000)(s+22.63+ jl4.55) 

(s+l .74)(s+1.712)(s+20)(s+22.63+ jl4.55)(s+22.63- jl4.55) ^ 

^ (s+22.63-il4.55) 

(S+80+ j52)Cs+80- j52)Cs+1660) ’ 

which after simplification is, 

s^(s+10)(s+1000) ,, 

( s+l. 712) (s+ 1660) (s+80+j52)CT+80-j52) ' 
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Using the graphical residue method on equation (6-2), 
the output due to a load disturbance applied as a unit step 
is, 



c = 


.990e”®°^sin(52t + 2.530) 


- .000965e 


-1.712t 




0.437e"^^^°^ . 




Pig. 


6-1 shows a plot of this 


response . 


The predicted 


sponse 


characteristics together with those 


obtained from 


g. 6-1 


are as follows: 










Predicted 


Ac tual 




Fall-time: < 


.01 sec 


.007 sec 




Overshoot : 


8.7% + 


8.6% 




Steady-state error: 


0 


0 




Oscillations : 




1 



The fall-time has been increased by 535% over that of 
the basic uncompensated system response. The overshoot has 
been decreased by 16.5% and no steady-state error results. 
This improvement in the regulation response did not result in 
any change in the desired original input-output response of 
the basic system. 

In order to check the accuracy of the correlation 
technique rules developed, the response of the representing 
function to a unit step input was determined using graphical 
residue methods, and was found to be, 

C (R p ) = .0488 + .709e"^°^sin(52t + 2.308) + 

.000745e"^*’^^^^ + .426e"^^^°* . 
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Fig. 6-1 shows a plot of this R. F. response. The difference 
in phase angles between the R. F. and the desired function 
is equal to (2.530-2.308) radians. Thus, the true value of r 
is equal to 0.222 radians or 12.7 degrees. The y obtained 
through the correlation rule, equation (5-8), was 13 degrees. 
The estimate of the R. F. was .079, while the actual value 
of was computed to be .0488. The accuracy of the 
correlation rules used in the representing function technique 
has been shown to be more than sufficient. The time saved 
in using this technique for determining the desired conditional 
feedback compensation was considerable. 
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7. Selection of disturbance input. 

The external load torque on the antenna due to wind 
gusting is a function of time with three independent 
variables, wind velocity, wind direction, and the tracking 
direction of the antenna. To evaluate the compensation, it 
is desirable to investigate its effectiveness in reducing 
the greatest expected disturbance error. Since the three 
variables affecting the disturbance output are independent, 
the maximum error will occur when the errors due to each 
variable are maximum. 

To investigate the effectiveness of the compensation, 

a wind velocity profile having the greatest expected severity 

was used. Prom consideration of data pertaining to wind 

2 

gust velocities , a representative velocity profile was 
developed which would produce the maximum expected disturbance 
error. This velocity profile is shown in Fig. 7-1. A steady 
25 miles per hour wind was assumed to exist prior to the 
gust. The wind velocity then increased smoothly to a peak 
of 60 miles per hour in one-half second and remained at this 
value for one-half second. The maximum gust velocity normally 
decreases at a slower rate than it increases during the gust 
build-up. Since only that part of the gust profile resulting 
in the largest disturbance error is of interest, the trailing 
edge of the actual gust profile need not be considered. 

^W. R. Gregg, Aeronautical Meteorology, pp. 269-274, 1930 
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Note: The wind velocity profile was developed as 

straight line segments approximating a 
smoothly increasing wind gust in order to 
utilize a function generator in simulating 
the wind velocity input to the analog 
computer simulation described in Appendix IV. 



Figure 7-1. Representative wind velocity profile. 
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The true wind direction was considered constant since 
wind direction changes are most likely to be negligible 
during the short interval of the gust being used. Since 
tiue wind has been considered constant during the one-half 
second of applied gusting, the change in relative wind 
direction is a function only of antenna rotation. The 
targets tracked by this antenna are earth-orbiting satellites 
which require relatively slow tracking rates. Thus, the 
change in relative wind due to antenna rotation is negli- 
gible during the one-half second gusting interval. During 
the gust, the antenna was considered to be in that position, 
relative to the true wind, which would cause the greatest 
disturbance torque for any given wind velocity. 

A plot of wind moment coefficient (C ) about the 

mg 

antenna gimbal center vs. angle of wind incidence has been 
prepared for the given antenna from wind tunnel tests made 
on a l/30th scale model. The data is reproduced in Fig. 7-2. 
The use of a ''spoiler,” fixed to the back side of the antenna 
around the periphery of the dish, was used to reduce maximum 
torque loading on the antenna by approximately 50 percent. 

The moment about the antenna gimbal center (M ) was 

s 

derived in Appendix III as 

Mg = 201.9 (ib-ft) , (7-1) 

3 

Guggenheim Aeronautical Laboratory, California 
Institute of Technology, Report 737, 1960. 
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B 



where C = dimensionless moment coefficient about gimbal 
center 

V = wind velocity (ft/sec), 
w ' ' 

For a maximum velocity of 60 MPH and a maximum C 

' mg 

of 0.189 (130*^ angle of wind incidence), the maximum M 

O 

computed from equation (7-1) is 137.3x10^ Ib-ft. This 

moment, M , is the external torque (T ) applied to the 
& ^ 

output azimuth drive shaft. 
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8. Evaluation of the conditional feedback compensation. 

The evaluation of the designed conditional feedback 
compensation is made on the basis of a comparison between 
the uncompensated disturbance-output response and the 
compensated system response. 

A step position disturbance input was used in Section 
6 to determine the design of the compensation since this 
input represented the most severe disturbance, though only 
theoretically possible. The compensation was designed to 
provide for an arbitrarily selected improvement in system 
response. These responses were shown in Fig. 6-1. 

A comparison of system response to a step velocity 
disturbance input was made to closer approximate the system 
response to wind gusting. These responses are shown in Pig. 
8 - 1 . 

In order to verify the degree of response improvement 
indicated as a result of the above two disturbance inputs, 
the system response to the wind gust profile of Fig. 7-1 
was obtained. The electronic analog computer simulation for 
this disturbance input was formulated in Appendix III. 

Appendix IV presents the analog simulation used to determine 
the system response to the wind gust. Both compensated and 
uncompensated system responses to the gust profile of Pig. 7-1 
are shown in Fig. 8-2a. 

In order to note the system response to a wind gust 
more severe than that derived in Section 7, and, in addition, 
one of less severity, wind gust profiles of one-half and 
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twice the time interval of the profile of Fig. 7-1 were 
applied. These responses are shown in Pig. 8-2b and Fig. 8-2c. 
The frequency response curves for the compensated system are 
shown in Fig. 8-3. 

Table I shows the disturbance response characteristics 
of the uncompensated and compensated systems for the 
disturbance inputs applied. The most significant parameter 
upon which to evaluate the compensated system is the maximum 
error caused by the wind gust. Table I shows this error, due 
to the wind gust of Pig. 7-1, to be 7.2 seconds of arc for 
the uncompensated system. With the addition of the conditional 
feedback compensation, this error was reduced to 0.028 seconds 
of arc. It is emphasized that this substantial reduction 
in antenna position error due to wind loading has been effected 
without change to the original input-output response. 
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DISTURBANCE RESPONSE CHARACTERISTICS 
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Phase Shift (Degrees) -78.0 -113 



9. Conclusions. 



This study has shown that there exists a compensation 
technique which can be successfully utilized to reduce 
antenna tracking errors due to wind disturbances without 
effecting a change in tracking input-output response. 

Section 2 introduced the conditional feedback compensa- 
tion method and has shown that this method permits design 
specifications on disturbance-output response to be met 
without change to the tracking input-output response. A 
specific conditional feedback compensation was developed to 
obtain a substantial degree of reduction of wind disturbance 
tracking error. The degree of disturbance error reduction 
is presented in Table I of Section 8. The representative 
wind velocity profile developed in Section 7 was applied as 
the disturbance input, and the resulting error decreased from 
7.2 seconds of arc without compensation to .028 seconds of 
arc with conditional feedback compensation. This represents 
a 257:1 reduction in disturbance error. In addition, this 
substantial error reduction was achieved with a faster response 
time. The response improvement of the compensated system is 
made more evident by noting that the other response parameters 
listed in Table I show substantial improvement over the uncom- 
pensated system. A significant disturbance error reduction 
might be expected when it is noted that bandwidth has been 
increased from 2 cps . to 98 cps. 

In order to determine the system response due to any 
given conditional feedback compensator, two root locus plots 
would have been required to place the roots of the closed loop 
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function. It was anticipated that several values of 
compensation would have to be investigated before a satis- 
factory compensator could be obtained. The representing 
function technique was originated for the purpose of expediting 
the above repetitive process. This technique was employed 
in the determination of the desired compensator. The represen- 
ting function technique proved to be accurate, expeditious, 
and simple to apply. A step-by-step procedure for its 
application is given at the conclusion of Section 5. 
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APPENDIX I 

COMPLEX ZERO COMPENSATOR 

4 

An electronic compensator, Fig. I-l, has been developed, 
utilizing an operational amplifier, which will essentially 
cancel the complex poles of a cascaded transfer function. 




Figure I-l. Complex zero compensator schematic. 



'‘l * ^2 * ^4 = ° 
«2 



(I-l) 

( 1 - 2 ) 



R, 



^2 ^3 



I3 - C2Vj^s 



- = ''l' K7 * ^ ^ c^s) 
eR2 

''1 " 



2 ^2 ^1^2^2® 



(1-3) 



G. A. Korn and T. M. Korn, Electronic Analog Computers, 
pp. 417-420, McGraw-Hill Co., 1956 
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1 . = 



E - V, 



R, 



= ^4 ^ ^5 



^5 = 



R 4 - ^^ 2 ^ R 3 ^ R 4 " ^3^^ 



ER, 



V^ = 



2 R^ + R4 + R3R4C2S 



•4 R3 + R4 + R3R4C3S 
Substituting ( 1 - 2 ), ( 1 - 3 ), and ( 1 - 4 ) into (I-l) gives 



( 1 - 4 ) 



C es + ? + 5 

R, + R^ + R,R^C^s R^ + R, + R^RX„s 



= 0 . 



^1 *'2 ‘' 1*'2 2 ‘ 



'3 *'4 ^'3 4 3 ' 



( 1 - 5 ) 

Letting R^ = R2 = R3 = R^ = 1 megohm, then equation ( 1 - 5 ) 
becomes 

e E 



Cj^es + 



2x10^ + C2sxl0^^ 2x10^ + C3sxl0^^ 



= 0 (1-6) 



e (2 x10^CjlS + Cj^C2s\lO^^ + 1) 
2xlO^ + C^soclO^^ 



-E 



2x10^ + C3sxl0^^ 



( 1 - 7 ) 



If C2 is made equal to C3, the denominators of ( 1 - 7 ) become 
equal, and ( 1 - 7 ) reduces to 



I = -C,C2Xl0^2(s2 * s . ) 



( 1 - 8 ) 



If the capacitors are valued in micro-farads, equation 
(1-8) becomes 



E _ ^ ^ /_2 . 2 

e ~ '^1^2 



- = -C,C^(s + ^ s + 



) 



1^2 



( 1 - 9 ) 
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Equation (1-9) is of the form of a constant times a complex 
zero, where and are picked so that the coefficients 
in this equation equal those of the complex poles it is to 
cancel . 

Depending on the values of capacitors used, it may be 
necessary to modify the circuit values slightly to suppress 
high frequency oscillations in the compensator amplifier. 
Reduction of the value of or placing a relatively small 
resistance (about 0.1 megolms) between and ground has 
been found effective. Such changes naturally reduce the 
cancellation effectiveness of the compensator, but to a 
degree generally within acceptable limits. 
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APPENDIX II 



GRAPHICAL RESIDUE TECHNIQUE FOR DETERMINING 
THE TRANSIENT RESPONSE EQUATION OP A SYSTEM 

Consider the disturbance-output transfer function 
of the uncompensated basic system: 

C _ 1 _ s^(s + 10) 

U 1 + F (s + 1.74)(s + 22.63 - J14.55) 
o 

^ (s + 22.63 + J14.55) * ClI-1) 

Let U(t) = A u(t), where A = 1. Therefore U(s) = — , and 

s 

equation (II-l) becomes 

P _ s( s -*■ 10) 

^ (s + 1.74)(s + 22.63 - jl4.55)Cs + 22.63 + J14.55) ' 

(II-2) 

The poles and zeros of equation (II-2) are plotted in Fig, II-l. 
The residue at (-1.74) is equal to 

TT ( vectors from the zeros to the pole ) 

TT ( vectors from the remaining poles to the pole) 

or, (r^ - - z^) r^^t 

[ri - ( " * ‘ 

Prom Fig. II-l equation (II-3) becomes 

(1.74 /l80 ° )(8.26 /qo ) 1.74 

e 

(25.45 /2 TT - 0^)(25.45 /(^^) 

-1 74 t 

or, -0.0222 e ^ ^ ^ (II-4) 



(H-3) 
/l8QO t 
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The residue at the complex pole ( - 22.63 + J14.55) 



is equal to: 



1 


^(- oc + jw^) - J 


i ! 


[(- « jw^) - Zj] 






E- « - ""i] 


1 


j^(- oc + jw^) - (- C 


X - jw^)] 



From Fig. II-2, equation (II-4) becomes 



( - oc + jw^) t 



(H-5) 



(26.9Zr)(19.25 (-22.63 * jl4.35)t 

e 

(25.4 5/TT-Q^ )(2 X 14.55 /Q^ ) 

(26.9 z^)(19.25 g (-22.63 + jl4.55)t 



(25.45/TT-0^)(14.55) 
-22.63t . 



2 j 



= 1.395 e 



sin(14.55t + vp + 0 ^ _ Tf + 0 ^) 



= 1.395 e ■^^'^^^sin(14.55t + 2.319) . (II- 6 ) 

Combining equations (II- 6 ) and (II-4) yields the complete 

transient response equation for the system ^ u(t) ; 

1 + F 



C = 1.395 e ■^^•^^*sin(14.55t + 2.319) - 0.0222 e 



By this same technique, the other disturbance-output 
transfer functions used in this thesis can be converted to 
their equivalent transient response equations for various 
input signals. A listing of these transient response 
equations for unit step and unit velocity inputs follows. 
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C _ 1 _ s^(s + 10) 

U “ 1 + ~ (s + 1.74)(s + 22.63 - jl4.55)(s + 22.63 + jl4.55) 



C(t) 

lu(t) 



1.395 e 22.63t sin(14.55t + 2.319) - 0.0222 e 



= 0.0520 e -22.63t sin(14.55t - 0.247) + 0.0127 e 



C _ 1 

U - 1 . 

= (s + 1.74)(s + 22.63 - il4.55)(s + 22.63 + il4.55)(s + 1000) 
(s + 1.412)(s + 80 - j52)(s +80 + j52)Cs +1660) 



C(t) 
lu( t) 



0.0488 + 0.709 e sin(52t + 2.308) + .000745 e 



+ 0.426 e 



•1660t 



= 0.00172 + 0.0488t + 0.00738 e sin(52t - 0.265) 

+ 0.000423 e _ q. 000257 e 



C 

U 



1 + F 



o J . 



1 " ^IS 



s (s + 10)(s + 1000) 

(s + 1.712)(s + 80 - j52)(s + 80 + j52)(s + 1660) 

= 0.990 e sin(52t + 2.530) - 0.000965 e "l’712t 

+ 0.437 e 



= 0.0104 e sin(52t - 0.0306) + 0.000564 e 



- 0.000263 e 



-1660t 
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APPENDIX III 

DEVELOPMENT OF DISTURBANCE TRANSFER FUNCTION 



The purpose of this Appendix is to develop the relationship 
between the disturbance input (U), a function of the disturbing 
wind velocity and the antenna position angle (9^) • The 

development of this transfer function was undertaken in two 
steps. First, the transfer equation of the disturbance torque 
(T ) with respect to V was determined. The transfer function 
of T to the antenna position angle (0.) was then determined 

0 i-/ 

and the two steps combined to give the transfer function of V 

w 

to 0L . 



Fig. 7-2 presents a plot of the non-dimensional moment 

coefficient (C ) of the antenna as a function of relative 
mg 

wind direction. Using Fig. 7-2 to determine C , the moment 

^ ^ mg’ 

about the antenna gimble center (M ) which is also the external 

b 

torque (T ), was found as follows:^ 



T=M=C qSr 
e g mg 



(III-l) 



where , 




the dynamic pressure 



(|| 2 ) 



1 



2 

S = 2830 ft , antenna dish frontal area, 

f = air density ( , 

f t4 

r = 30 ft, antenna dish radius. 

The data of Fig. 7-2 was taken at a temperature 

and a barometric pressure of 742 mm Hg . 

The air density ( P ) is defined as, f = — ^ 

gKl 

^J. H. Dwinell, Principles of Aerodynamics, pp. 
51-52, McGraw-Hill Co., 1949 



of 20°C. 



(III-2) 

25-26, 
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where , 



p = absolute pressure ( 

ft^ 

g = 32.2 2 , acceleration due to gravity 

sec^ 

R = 53.3 — , gas constant 
°R 

T = absolute temperature (°R) . 

The temperature of 28°C is equal to an absolute temperature 
of 541. 9°R. A pressure of 742 mm Hg. is equal to 2075 Ibs/ft 
at 28®C . Substituting these constants into equation (III-2) 
gives , 

f = .002235 Ih-sec^/ft^ . 

Solving for q in terms of and substituting into 
equation (III-l) gives, 

T = 201.9 C V ^ ft-lbs . (III-3) 

e mg w 

Equation (3-18) of Section 3 shows the relationship 
between the error signal (X) and the external load torque 
(T^) with respect to the position output angle (0^^) of the 
hydraulic motor. Fig. 3-2 shows the hydraulic system configu- 
ration. Thus, the actual input to the hydraulic function (Gj^) 
of Pig. 3-2 for an applied external wind torque is the sum 
of the error signal (X) and some function of external torque. 
The hydraulic motor output angle (0^^) can then be written as, 

0^ = X + T^ (III-4) 

where G^^ is some transfer function which carries external 
torque into an error position to be applied to the hydraulic 
transfer function (G^^) . Letting X in equation (III-4) be equal 
to zero, then 0^ = G^G^^ T^ . (III-5) 
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I 





Substituting X equal to zero into equation (3-18) gives, 



0 , 



(-) Kj_ (L.^s) 



dm N(Jj^s +fj^s+Kj^) 



dm s+(L+^s)J s^+(L+^s)K, (J, s^+f.s) 

D m D L ^ L L 

“din "^N2(j^s2+fj^s+Kj^) 



(III-6) 



Substituting equation (III-6) into (III-5) gives, 






(-) (L+|s) 



s []dm^ N(Jj^s^+fj^s+Kj^) + (L+^s)NJ^s(Jj^s^+fLS+Kj^) + (L+|s)_I^( Jj^s+fj^)~| 

N 



(in-7) 



Making the same assumptions as in Section 3, that is, K. is 

fL 

very large and — ^ = 0, equation (III-7) simplifies to. 



^iS = 



N'-CJ +J, ) 
m L 

n2 



(-) (L+|s)NB/V(N^J^+Jj^) 



s 


s^+BLs+B dm^N^ 






V V(n2j +j^ ) 

— m L _ 





(III-8) 



Dividing equation (III-8) by as given in equation 

gives, „ 

. (GlGh) -(L.|s) 

G. ■ NK dm 
h V 

Fig. III-l shows the final block diagram of the 
system with the added complex zero compensation (G^) 
tachometer feedback (G^-) described in Section 4. In 
the application of through the synthetic transfer 

is represented. 



(3-21) 

(III-9) 

hydraulic 

and 

addi tion , 
f unc tion 
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Figure III-l. Block diagram of final hydraulic system. 



Assuming that 0^^ is due to the external torque alone, 
( ie . X ' = O) , then 




T 

e 









(III-IO) 



where, = K^s, the tachometer feedback 

G = K X (complex roots of G. ), the complex zero 
o o h ’ 

compensator added in Section 4. 

Substituting equation (III-8) into (III-IO) and simplifying 
gives. 



©L _ (s+BL/V) / 



T 

e 


s 


s2+BL/Vs+ B dm2 ^2 




1 + KfKoKyBN^ dm 






V(N^J +J, ) 

_ m L _ 




V(J,+N^J ) 



(III-ll) 

Substituting the hydraulic flow coefficients listed in 
Section 3 and the gain constants developed in Section 4 into 
equation (III-ll) gives. 



_ 1.51x10 ^(s->-42) rad/lb-in . (III-12) 

^e s(s^+42s+3228) 



77 



Combining equations (III-12) and (III-3) enables the 



solution of the output position angle of the antenna for 
any given wind velocity. Thus, for a constant , 



_ 201.9 Cmy x 1 . 51xl0~^( s+42) 
Vw^ s(s^ + 42s + 3228) 



(III-13) 
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APPENDIX IV 



ANALOG SIMULATION OF CONDITIONAL 
FEEDBACK COMPENSATED SYSTEM 

Fig. 5-2 shows the block diagram of the complete 
compensated system under study. It has been shown that only 
the response of the system to a disturbance input is of 
interest since the response to the tracking input (R) is 
knoim to remain unchanged with the addition of the conditional 
feedback compensation. Fig. IV-1 shows the block diagram of 
the compensated system with the tracking input equal to zero. 




Figure IV-1. System with conditional feedback and R = 0. 

The transfer blocks of Fig. IV-1 were defined as follows 

in Section 4: 

F = 37(s + 1.7)(s + 20) 

° s^ (s + 10) 

P _ 21.06(s + 96)(s + 200) 

^2 (s + 20)(s + 1000) ’ 

and where U is the disturbance input, antenna 

position due to an applied wind velocity. Equation (III-13) 

has shown the transfer function from wind velocity (V ) to 

' w 
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0j^ to be , 



U 



V 



w 



01 

V 



201.9 C 



mg 



X 1.51 X 10 ^(s + 42) 



w 



s(s + 42s + 3228) 



The analog simulation of the compensated system of 
Fig. IV-1 was undertaken in order to determine the effects 
of the added conditional feedback compensation on the system 
response to a simulated wind velocity disturbance profile. 

The following is a list of equipment utilized in the 
simulation of the compensated system: 

3- Conner analog computers, Mod. 3000, 10 channels 
per computer. 

1- Brush dual-channel D.C. amplifier. 

1- Brush dual-channel pen recorder. Mod. BL-202. 

1- Conner variable-base function generator. Mod. 
3750. 

1- Conner electronic function multiplier. Mod. 3730. 

1- Conner null, voltmeter. Mod. 5002. 

1- Hewlett Packard low frequency function generator. 
Mod. 202A. 

1- Moseley Autograph x-y recorder. Mod. 2-S. 

The simulation of the disturbance input function was 
undertaken in two parts. First, the disturbance torque (T^ 
was simulated using equation (III-3), 



T = 201.9 C V ^ ft-lbs 
e mg w 



The representative wind velocity profile of Fig. 7-1 
was simulated as a function of time using a variable-base 
function generator. The simulated profile was generated using 
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I 

( 

!<! 



i 



a step velocity input as a time base and an eleven straight- 
line segment approximation of the wind velocity profile. The 
constant initial velocity of 25 MPH was then added to the 
output of the function generator. The resultant was the 
simulated profile of Fig. 7-1. This velocity profile was 
squared using the function multiplier. Since was considered 

constant during the interval of the gust, the squared velocity 
was multiplied by the appropriate constant to give T^ . 

Equation (III-12) gives the transfer function between 
T and the output position due to wind loading (0^). The 
analog simulation of this transfer function requires an 
integrating output stage. This equation represents an open- 
loop transfer function. Therefore, there is no feedback loop 
around this output stage. Since the input to this integrator 
is non-zero, this integrating stage would saturate. The 
simulation was modified to include this integrating stage 
within the feedback loop of the servo system, resulting 
in an output from the disturbance simulator of the velocity 
(0j^) . Modifying equation (III-12) gives, 

y_ _ ^ (s * 42) 

e e (s^ + 42s + 3228) 

- 8 

where is equal to 1.51x10 

The velocity 0j^ was added to the velocity output of the 
transfer function and integrated to obtain the system 
position output (C) . 

Fig. IV-2 shows the complete schematic of the analog 
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computer simulation used. Values of components used as well 
as the primary outputs are indicated on the figure. 
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